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Abstract 


n .™V? dC devel °P m « nt and application program for the Langley Aerothermodynamic 
Upwind Relaxation Algorithm (LAURA), with emphasis directed toward support of the 
, 4 e rr, FUght Ex P eriment CAFE) in the near term and Aeroassisted Space Transfer Vehicle 
(ASTV) design in the long term is reviewed. LAURA is an upwind-biased, point-implicit 
relaxation algorithm for obtaining the numerical solution to the governing equations for three- 
dimensional. viscous, hypersonic flows in chemical and thermal nonequilibrium. The algorithm 
is derived using a finite-volume formulation in which the inviscid components of flux across 
cell walls are described with Roe’s averaging and Harten’s entropy fix with second-order correc- 
tions based on Yee s Symmetric Total Variation Diminishing scheme. Because of the point- 
implicit relaxation strategy the algorithm remains stable at large Courant numbers without the 
necessity of solving large, block tri-diagonal systems. A single relaxation step depends only on 
information from nearest neighbors. Predictions for pressure distributions, surface heating, and 
aerodynamic coefficients compare well with experimental data for Mach 10 flow over an AFE 
wmd tunne model. Predictions for the hypersonic flow of air in chemical and thermal none- 

nhrJTrt 111 (ve oc ? ty ~ 8917 m/s ’ aitltude = 78 km.) over the full scale AFE configuration 
obtained on a multi-domain grid are discussed. * 
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SIMULATION CAPABILITIES 

PROGRAM LAURA (LANGLEY AEROTHERMODYNAMIC 


05 ^ 

O CD 


co -3 

>- z 



-c* 

0 ) 


CD 

'0 



CD 


05 

CL 


— 1 

CO 


>* 

k. 

1 — 


05 

1 — 


ID 

" — ' 


C 

=3 

£ 


O 



CD 

*r~ 



.O 

“O 



05 

c 


CD 
-♦— * 

0 

cr 

CO 


CD 

a 

_Q 

C 

Z> 


O 

L. 

H 

Z 

O 


To 


05 

0 

ID 

c 

E 

03 

0 ) 

E 

05 

CD 

sz 

CO 

— J 

0 


CD 

CL 

E 

CD 

o oo 
5 CNJ 


3 < 

~ CO 
J Q < 
= Z 

cr c 

CD — 

§ -o 

Z S 
c 

- 7 =: CD 

i I 

Q) O 


CO > 

y « 

DC 

LLl C 


486 


Multi-Domain Capability, Structured Grids 

Ideally Suited For Parallel, Asynchronous Iteration Using Unstructured Grids 



COMPUTATIONAL GRID ON SURFACE AND PLANE 
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PRESSURE AND HEATING DISTRIBUTION FOR AFE 

Mach 10 a = 5° 
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PRESSURE AND HEATING DISTRIBUTION FOR AFE 

Mach 10 a = 0° 
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S/L 



PRESSURE AND HEATING DISTRIBUTION FOR AFE 

Mach 10 ot = -5° 




DETAIL OF AXIS AND NO AXIS GRIDS 




HEATING DISTRIBUTIONS FOR AFE STING 

Mach 10 a = -5° 
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PRESSURE CONTOURS IN THE PLANE OF 
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MACH NUMBER CONTOURS IN THE PLANE OF 



495 


COMPUTATIONAL GRID ON SURFACE AND 
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TRANSLATIONAL TEMPERATURE CONTOURS 
IN THE PLANE OF SYMMETRY 



VIBRATIONAL TEMPERATURE CONTOURS 
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ELECTRON NUMBER DENSITY CONTOURS 
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TEMPERATURE DISTRIBUTION ACROSS THE 

SHOCK LAYER 




